Renormalization of quark axial current in the chiral potential model 
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Non-conserved composite operators like the quark axial current have divergent matrix elements 
therefore must be renormalized. We explore how this can be done in quark model calculations where 
the systematic procedure of dimensional regularization and minimal subtraction is not applicable. 
We propose a most natural and convenient regularization scheme of cutting the intermediate quark 
states over which we sum in loop diagram calculations at a certain energy. We show that this scheme 
works perfectly for the quark axial current and we obtain the quark spin contribution to the proton 
spin: Au = 0.82, = —0.43, As = —0.10, which is in excellent agreement with experiments. 

PACS numbers: 12.39.Ki, 12.39.Fe, 12.39.Pn, 13.88.+e 



The expression nucleon spin "crises" denotes the find- 
ings of the European Muon Collaboration (EMC) [Q that 
a small proportion of the nucleon spin is carried by the 
quark spin and that strange quark polarizes significantly 
in the nucleon. This has been under hot debate for over 
ten years but one obtained not yet a fully satisfactory de- 

', scription (for a review of the nucleon spin problem, see, 
e.g., lH). It should be emphasized that the "crises" is 

\ not for the fundamental theory of quantum chromody- 
namics (QCD), in the viewpoint of which the nucleon is 
a complicated object of quarks and gluons and quarks 
do not necessarily carry most of the nucleon spin. The 
"crises" is, however, for the naive SU(6) quark model 
which is quite successful in many aspects but neverthe- 
less attributes all the nucleon spin to constituent quarks. 

\ To explore whether this "crises" is real, i.e., whether the 
SU(6) model could be taken as a good lowest order ap- 
proximation for the nucleon, a natural way is to start 
from the SU(6) wavefunction and study whether we can 
explain the experimental result of the nucleon spin con- 
tent by going to higher orders. 

In the past years there has been countless work us- 

, ing quark models along this direction (for references see 
1^), but the problem is not really resolved. The obsta- 
cle is that the quark axial current, which is the operator 
for defining the quark spin contribution to nucleon spin, 
is a non-conserved composite operator. Therefore when 
we go to higher orders divergent matrix element will be 
encountered and the quark axial current must be renor- 
malized. But unfortunately, the usual renormalization 
schemes are not applicable in quark model calculations, 
where instead of divergent loop integrations over the con- 
tinuous momentum, we encounter divergent summations 
over the discrete quark excited states whose wavefunc- 
tions are obtain numerically. To explore how one can 
renormalize a composite operator in quark models and 
what would be the result of the renormalized quark axial 
charge (i.e., the quark spin contribution to nucleon spin) 
are the aims of this paper. 

In the following wc will first construct the bare matrix 



element of the quark axial current and demonstrate its 
divergence, then we explore how we can renormalize the 
bare quantity by systematically subtracting a divergent 
part and obtain a (finite) physical result. 

The quark spin contribution to the nucleon spin is de- 
fined as the quark axial charge of the nucleon: 

{Ps\tpql^l^i^q\ps) = Upsl^-i^Up, ■ A,, (1) 

where q — u,d,s and Ups is the nucleon spinor. An equiv- 
alent but more suitable expression for model calculation 
is 

Here \p+) is a nucleon state with positive momentum and 
polarization along the third direction. We are going to 
adopt Eq. (H) to pursue a perturbative calculation of 
in a chiral potential model. Our model Lagrangian is 

£ = 5'(r) - 7V(r)]V' - 

^i^[S{r){a + Z7'AV.) + (a + i-i^X'<j,,)S{r)]i^ + 

\{d,af^\{d,<^.f-\mla^-\mjc^l (3) 

The model Lagrangian is derived from the a model in 
which meson fields are introduced to restore chiral sym- 
metry 1^]. The flavor and color indices for the quark 
field "0 a-re suppressed; the scalar term S{r) = cr + m 
represents the linear scalar confinement potential cr and 
the quark mass matrix m; V{r) — —a/r is the Coulomb 
type vector potential and F^=93MeV is the pion decay 
constant, a and 1/)^ (i runs from 1 to 8) are the scalar 
and pseudoscalar meson fields, respectively and A,; are 
the Gell-Mann matrices. The quark-meson interaction 
term of Eq.(^ is symmetrized since the mass matrix m 
does not commute with all for different quark masses. 

At zeroth order the nucleon is taken as the usual SU(6) 
three-quark ground state of the Hamiltonian 
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H,^ = / Sxil^^[d ■ -d + l3S{r) + V{r)]ij. 



(4) 



The diagrams for the numerator and denominator of Eq. 

up to second order are shown in Figs. 1 and 2 respec- 
tively. 
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FIG. 1. Feynman diagrams for the matrix element 
{N\ ^ dPxTl)^'^^^il)\N) up to second order; a cross on the 
quark line denotes the quark axial vertex; A is of the ze- 
roth order, B is the renormalization counter term, C and D 
are vertex and exchange diagrams respectively. The meson 
line in C can be a tt, r;, or cr (while the intermediate quark 
is u or d), or & K (while the intermediate quark is s); the 
meson line in D can only be a tt, 7; or a. 
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FIG. 2. Feynman diagrams for the normalization {N\N) 
up to second order; A is of the zeroth order which is simply 
unity, B is the meson exchange diagram. The meson line in 
B is a TT, rj, or a . 
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FIG. 3. Charge renormalization condition; a cross on the 
quark or meson line denotes the zeroth component of the 
vector current vertex. 

By computing Figs. 3C and 3D we can determine the 
renormalization constant Z2 , which is then to be used in 
Fig. IB. 

Next we remark that at zero momentum transfer, the 
exchange diagram Fig. ID is actually the product of Figs. 
lA and 2B, therefore the sum of Figs. lA and ID over the 
normalization Figs. 2A and 2B is just Fig. lA. (However 
this is not true when we calculate the axial form factor 
at finite momentum transfer.) Thus only Fig. IC is left 
to be evaluated together with Fig. 3C and 3D. 

The essential ingredients needed for calculating these 
diagrams are the quark and meson propagators. The me- 
son propagator given by the Lagrangian of Eq. (||) is the 
free propagator: 



{2n) 



(7) 



Since the non-perturbative confinement is included in Hq 
the quark propagator has to be obtained numerically, and 
in practise we have to work with time-ordered perturba- 
tion theory. We write the solution of Hg as 



\ua{x)aa 



^vp{x)hl, 



(8) 



We first discuss how to determine the renormalization 
constant Z2. The mass-shell renormalization scheme is 
not applicable here, since our unperturbed quark basis 
are confined wavefunctions. But we can still use the 
charge renormalization condition. The conserved elec- 
tromagnetic current of the Lagrangian of Eq. (^) is 
f = Y^qjq + 34,, where and are the quark and 
meson current respectively: 

fq = Qqi^^r^g, 

= e(0ia^02 - 029^ + cbidf^cb^ ~ cb^df^M- (5) 

The charge renormalization condition is to require that 
for a quark state 



where Ua{x) = e ''^Ua{x)Ta, vp{x) 



'^^Vp{x)Tp] T 



{q\ / d^xf{x)\q)^Qq. 



(6) 



Up to second order this is shown in Fig. 3. 



is the flavor wavefunction and Ua{x) and Va{x) are the 
spatial wavefunctions. The quark propagator is then 

D{xi,X2) = 0{tl - t2)'^Ua{xi)Ua{x2) - 

a 

Vl3{xi)vp{x2)- (9) 

13 

It can be shown by carrying out the time and energy 
integration that apart from an isospin factor, Figs. 3C 
and 3D yield the same expression. Therefore we define 
for Figs. 3C and 3D the pure space-time amplitudes: 

^0 = J d^xd'^Xid^X2A{x2,Xi) x 

Uf{x2)T^D{x2,x)j°D{x, xi)T^Ut{xi), (10) 

where Ui and Uf are the initial and final quark state 
respectively, the vertex function Ttt^k,?] = S{r)j^ and 
F(j = ~iS{r). All flavor wavefunctions are dropped out 
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and it should be understood that for cf) = K the inter- 
mediate quark is the s quark and otherwise is the u or d 
quark. Accordingly for Fig. IC we define 

Uf{x2)T^D{x2,x)^^^^D{x, xi)r^Ui{xi). (11) 

Now we can express the renormalization constant Z2 and 
the axial charge Aq in terms of and multiplied by 
spin and isospin factors which are calculated straightfor- 
wardly: 

= 1 - (35, + 2Bk + + B,), 

Zi = l-iBK + ^B^ + B,), (12) 

4 2 4 4 

A„ = -JrZ^ + -A^ + -Ar, + -A, 

Ad - -IIrZ^ + - - ^A, 

A, = 2Ak. (13) 

Another useful relation is for the nucleon axial charge 

gA - A„ - Arf - ^fnZ^ - + ^A, + ^A, (14) 

In Eqs. ( p^ -(p^) we have assumed equal masses for 
u,d and for Tr",?!^, therefore Z2 = Z2 which is just a 
statement of SU(2) symmetry (If SU(3) symmetry is un- 
broken Z2 would also be the same as Z^' ); 5/3 fa is the 
zeroth order values of gA and fu is a relativistic reduction 
factor. 

Since A^ and B^ correspond to loop diagrams, they 
would naturally be divergent. If the quark axial current 
were conserved, then the divergence of A^ and B^ would 
automatically cancel in Eq. ( p^ ) and we get a finite 
result for A^. This is the phenomenon that conserved 
operators do not need extra renormalization besides the 
usual renormalization of mass, charge and wavefunction. 
However the quark axial current is not conserved, the di- 
vergences of Afp and B^ will not cancel and the naively 
obtained in Eq. (Oh is divergent, which is just the 
general case that a composite operator has divergent ma- 
trix element and needs extra renormalization, i.e., we are 
going to subtract a divergent piece from A^ and si- 
multaneously and leave a finite part in Eq.(|l3|). This fi- 
nite leftover would depend on how we regularize A^ and 
B^ and how much is to be subtracted as the divergent 
part. This is the renormalization scheme dependence. 

In the usual plane-wave perturbation theories, we have 
a most powerful and systematic renormalization scheme 
of dimensional regularization and minimal subtraction 
(MS) or modified minimal subtraction (MS), but it is ev- 
idently not applicable here. We must find a proper renor- 
malization scheme for the quark model calculations. The 
divergent integration over the (continuous) momentum 



in the plane- wave perturbation theories is here contained 
in the summation over the (discrete) quark intermedi- 
ate states. Thus we propose to regularize by cutting the 
summation at a certain energy, which is analogous to the 
lattice regularization by a finite lattice spacing. In the 
following we explore how this scheme works in practice. 

Our model parameters are listed in Table I. As we 
demonstrated in a recent paper , A^ is actually very in- 
sensitive to the model parameters in the above described 
energy-cutoff regularization scheme. This property is 
also found for B^. Therefore we do not take too much 
effort in choosing the parameters except for a fit to the 
nucleon and A masses with meson exchange potentials. 

Figs. 4 and 5 give the numerical results of B^ and A^ 
as a function of the maximum energy up to which one 
sums the intermediate quark states. The divergences are 
clearly seen. The problem is now how to determine the 
cutoff. In general the cutoff might vary from one dia- 
gram to another, without any guidance we would be at 
lost. Fortunately, besides Ag — 2Ak, Eq. ( p^ provides 
another clean relation: 

A„ + 4Arf = 10A,. (15) 

Thus with the experimental results: 

A, =0.80(6), Arf = -0.46(6), A, - -0.12(4), [7] 
A„= 0.82(6), Arf = -0.44(6), A, = -0.10(4), | (16) 

we can determine the cutoff for Ak and Aj^ using Fig. 
5. It is amazing to notice that the cutoffs needed for 
Ax and A^^ are roughly the same, i.e., independent of 
whether the intermediate quark is u,d (for A,) or s (for 
Ak) and also independent of the meson masses. This re- 
minds us that in the MS or MS scheme, the subtraction 
is independent of the mass parameters. We will therefore 
choose a ^-independent cutoff. It is also very interesting 
to notice that the cutoff value needed here roughly equals 
the inverse of the lattice spacing in lattice QCD cal- 
culation of Aq. To reduce the number of parameters, we 
choose the cutoff for to be the same as a^^ = 1.74GeV 
in j^]. But there remains still one important question to 
be asked: should the cutoff for and be the same? 

TABLE I. Model parameters and basic model predictions; 
the units for mass, a, and c are MeV, MeV-fm, and MeV/fm, 
respectively; center of mass corrections are made for the quark 
core contribution using the Peierls-Thouless method [5]. 
is the zeroth order nucleon axial charge, the full qa is obtained 
later. 
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It might be taken for granted that they are the same, 
(the simplest example is that we use the charge renormal- 
ization condition to determine the renormalization con- 
stant and then calculate the matrix element of the charge 
operator itself). However, we call special attention here 
that this is not necessarily the case. The problem is the 
lack of Lorentz covariance in the model. 

It is actually not difhcult to understand this point: If 
a theory does not respect Lorentz covariance, then the 
renormalization constant determined with the time and 
spatial component of the electromagnetic currents are 
possibly different. Therefore if we calculate the matrix el- 
ement of the charge operator but use the renormalization 
constant determined through the renormalization condi- 
tion for the spatial component of the vector current, then 
we might not obtain the physical charge (in some cases 
the result may still be finite, but this is not enough). 
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FIG. 4. Plot of as a function of the maximal energy 
up to which the intermediate states are summed over. 



stant determined for the charge might not be suitable for 
calculating the axial charge. But since we have presently 
no better method, we must allow different cutoffs for 
and B(j,. 

After describing the formalism, we are now in the po- 
sition to see how this renormalization scheme can ex- 
plain the experimental results. As we explained the cut- 
off for ^0 is chosen as the inverse of the lattice spacing 
a-^ ^ 1.74GeV §. Using Eq. (|l|) and Figs. 4, 5 and 
requiring gA = 1.257, the cutoff for is determined to 
be 0.734GeV. Then combining Figs. 4, 5 and using Eq. 
(|l3|), we finally obtain: 

A„ = 0.823, Arf = -0.432, A, = -0.104. (17) 

By adjusting only one free parameter (the cutoff for 
Bcf,) to fit gA, we reproduced the experimental results of 
^u,d,s perfectly. This can be regarded as a success of 
our model and renormalization scheme. Since we have 
started from the SU(6) wavefunction as the zcroth order 
approximation, we could say that the spin "crisis" for the 
naive SU(6) model is not real. 

We end our discussions by emphasizing three points: 
(1) Non-conserved operators have divergent matrix ele- 
ments and even in model calculations they must be con- 
sistently renormalized. (2) A natural and convenient 
renormalization scheme in quark model calculations is 
to cut the summation over the intermediate quark states 
at a certain energy; a fit to the experimental data reveals 
that this approach shares the same advantage of mass- 
parameter-independence as in the MS or MS scheme. (3) 
However, due to the lack of Lorentz covariance in quark 
models, the cutoffs for operators of different Lorentz type 
are not necessarily the same. 
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FIG. 5. Plot of A,f, as a function of the maximal energy 
up to which the intermediate states are summed. 

In the present model, Lorentz covariance is violated by 
the static potentials, therefore the renormalization con- 
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